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ABSTRACT 

Context. In the interpretation of stellar jet observations, the physical parameters are usually determined from emission line ratios, 
obtained from spectroscopic observations or using the information contained in narrow band images. The basic hypothesis in the 
interpretation of the observations is that the emitting region is homogeneous along the line of sight. Actually, stellar jets are in general 
not homogeneous, and therefore line of sight convolution effects may lead to the main uncertainty in the determination of the physical 
parameters. 

Aims. This paper is aimed at showing the systematic errors introduced when assuming an homogeneous medium, and studying the 
effect of an inhomogeneous medium on plasma diagnostics for the case of a stellar jet. In addition, we explore how to reconstruct the 
volumetric physical parameters of the jet (i. e., with dependence both across and along the line of sight). 

Methods. We use standard techniques to determine the physical parameters, i. e., the electron density, temperature and hydrogen 
ionisation fraction across the jet, and a multi-Gaussian method to invert the Abel transform and determine the reconstructed physical 
structure. 

Results. When assuming an homogeneous medium the physical parameters, integrated along the line of sight, do not represent the 
average of the true values, and do not have a clear physical interpretation. We show that when some information is available on the 
emissivity profile across the jet, it is then possible to obtain appropriate derivations of the electron density, temperature and ionisation 
fraction. 

Key words. ISM: Herbig-Haro objects - ISM: jets and outflows - Techniques: image processing - Methods: data analysis - Stars: 
pre-main sequence - winds, outflows 



1. Introduction 

Stellar outflows are generally observed as chains of bright knots 
with a characteristic spectrum that mainly includes forbidden 
emission lines (see, e. g., the review by Reipurth & Baliy |200TI) . 

Different approaches are used to extract the information on 
the outflow excitation conditions from the observed emission 
lines. To calculate the electron density n e the standard method 
involves the line ratio [SII]/16716//16731, which is strongly de- 
pendent on n e and less sensitive to the electron temperature 
T (see, e. g. Osterbrock 1 1989b . The [SII] ratio, together with 
the ratios between lines with different excitation temperatures 
(e. g. [OI]/15577//16300) is used to determine both n e and T. 
Moreover, the [SII], [OI]X16300, 6363 and [NII]X*6548, 6583 
emission lines are often used to calculate, apart from n e and 
T, the hydrogen ionisation fraction Xh and the total density 
n H (=n e /x H ) (Bacciotti & Eisloffel [T9991 hereafter BE99). 
Additionally, several line ratios in the optical and near infrared 
wavelength range (e. g. Brugel et al. 119811 Pesenti et al. 2003; 
Nisini et al. 120051 Podio et al. 120061) can be combined, with 
weights representing their respective errors (Hartigan & Morse 
120071 hereafter HM07), to determine the physical parameters of 
the HH objects. 
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Another possibility consists of obtaining predictions of emis- 
sion line intensities from plane parallel shock models (e. g. 
Hartigan et al. 119871 Hartigan et al. 119941 Lavalley-Fouquet et 
al. 1200 0, Pesenti et al. 2003) or from numerical simulations that 
include a detailed treatment of the ionisation and recombination 
evolution of several different species (e. g. Raga 2007 ), and com- 
pare the results with observations. 

It is important to develop strategies to properly extract the 
maximum possible amount of information from the observed 
line ratios. Convolution, instrumental and observing conditions 
(i. e. seeing) affect the derivation of physical parameters from 
observations. Projection or convolution effects are an obvious 
consequence of the lack of information on the spatial distribu- 
tion along the line of sight. Additionally, if the cooling region is 
not resolved the observed emission lines come from regions with 
different excitation conditions along the main jet axis. Actually, 
also in observations where the cooling region is well resolved 
the determined parameters are a convolution of the volumetric 
parameters along the line of sight, i. e. across the jet. 

Only for an homogeneous medium the mean values of n e , T, 
xh (integrated within the beam size and along the line of sight) 
will coincide with the corresponding values calculated from the 
observed line ratios. In the case of an inhomogeneous medium 
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it is necessary to understand what relation holds between the 
derived and the volumetric mean physical parameters. 

This problem has been studied in detail in the context of 
stellar atmospheres by several authors (e. g. Doschek 1984, 
Alml eakye t al. [1989l hereaft er AB S 89, Brown et al. [T99"Tl Judge 
et al. IT9971 Mcintosh et al. [19981 Judge l2000l As they show 
the information contained in a set of line ratios may be used to 
determine the function £{n e ), defined as the "emission measure 
differential in density", given by nl(z)dz = £(n e )dn e (e. g. Judge 
2000), where z is the position across the stellar atmosphere. If 
several line ratios are available, and if the geometry of the sys- 
tem is known, it becomes possible to infer the dependence of the 
density n e along the line of sight. 

Despite its importance, this problem has not been studied 
in similar detail in the context of stellar outflows. As far as we 
are aware, only Safier (1992) studied the influence of inhomo- 
geneities on the density diagnostics in outflows from T Tauri 
stars, applying a simple model of an isothermal medium with a 
power law radial dependence for the electron and total densities, 
and showing that different ratios may be used to obtain informa- 
tion on the density structure. 

Recently, observations with information on emission profiles 
across the jet have been presented by several authors (e. g. Beck 
et al. 120071 Bacciotti et al. |2000l HM07, Coffey et al. |2U07l ). If 
the emission line profiles are limited to a few observational val- 
ues across the jet we will show that this information can be used 
to derive the values of the physical parameters n e , T, and xh as- 
suming a Gaussian profile for the emission lines. If the emission 
profiles across the jet are well sampled, we will show in this pa- 
per how standard tomographic techniques (e. g. Craig & Brown 
1986 Brown 1995) may be applied to reconstruct the volume 
emission line intensities (i. e., the emission coefficients) from the 
values integrated along lines of sight, which can then be used to 
determine the « e , T and xh cross sections of the jet. 

The paper is organised as follows. In Section 2, the effect of 
stratification on the interpretation of the electron density, temper- 
ature and hydrogen ionisation diagnostics is discussed. Section 
3 presents a simple technique for obtaining information on the 
radial structure of the jet, and an example of its application to 
data of the HH30 jet. Finally, Section 4 summarises the obtained 
results. 



2. Convolution effects on jets 

2. 1 . Convolution effects for a non-resolved jet cross section 

We consider an axisymmetric jet in a three - dimensional carte- 
sian reference frame. The z-axis of this system is taken along 
with the main axis of the jet. The observer is looking at the jet 
along the y-axis, and the xz -plane coincides with the plane of the 
sky (see Fig.[T|). The cylindrical radius r is the distance between 
an arbitrary point and the z-axis. 

The jet is observed in the [SII]X16716, 6731, 
[OI]^6300,6363 and [Nil] ^6548, 6583 emission lines. 
These lines give three ratios which depend on the three variables 
to be determined: n e , T, x#. These emission line intensities are 
observed at a certain number of points x,- (corresponding to the 
pixel size) in the plane of the sky. 

In the following, capital letters will be used to design quan- 
tities integrated and convolved (according to the nature of the 
variable) along the line of sight, and lower case letters for vol- 
ume quantities. For instance, i represents the volume emissivity 
(defined by eq. IA.ll ). while / is the intensity (i. e. the emissivity 
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Fig. 1. Schematic representation of a cross section of a stellar jet. 
The jet is moving in the direction perpendicular to the plane of 
the image (the z-axis, not shown in the Figure), and is assumed 
to be axisymmetric. The observer is looking at the jet along the 
y-direction, and the observations are represented by a series of 
points along the x-axis, each one corresponding to the pixel size 
and obtained integrating the volume intensities along the line of 
sight and on the beaming size AxAz. 



n e represents the electron density across the jet, while N e is the 
electron density obtained convolving n e on a volume element. 
Only the electron temperature is designed by T in both cases. 

To clarify the effect of inhomogeneities on plasma diagnos- 
tics, first some results from ABS89 are reviewed. The ratio Gn 
between two observed line luminosities is given by: 



G\2 = 



JhdV 
j i 2 dV' 



(1) 



where dV is the volume element. Volume and integrated ratios 
have the same value if the medium is homogeneous. 

In general the line ratios also depend on T and xh, but the fol- 
lowing analysis is focused on ratios between lines emitted from 
the same species (independent of xh) and nearly independent of 
T (e. g. the red [SII] ratio). Assuming a medium with a uni- 
form electron density « e .o, eq.[T]leads to Gn = hlh - gi2(»e,o)- 
Inverting this equation it is possible to determine the density as 

n e ,o = gniGn). 

If the medium is inhomogeneous, the previous formula may 
also be applied to give a "spectroscopic electron mean density" 
defined as (ABS89): 



< n e >= 



JhdV 



J i 2 dV 



(2) 



The volume emissivities for the line j may be written as (e. g. 
ABS89. Osterbrock 1989 Safier [19921) : 



(3) 



1 J l+ njnj ' 
This expression is exact for 



a two-level atom. In this case, 
is the Einstein-A coefficient 



n j - A-2,\l{Ci,\ + C\y) is the "critical density" and ^ 

4^V2,i(Ci,2/xH)(n s /n H ), where A 2; i 
for the transition, C\ 2 and C2,i are the collisional excitation and 



integrated along the line of sight, with units: erg s cm sr ); de-excitation rates respectively, /iv 2j i is the transition energy, and 
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n s /«H is the population fraction of the species s. A fit with eq.[3] 
is also a very good approximation for many forbidden emission 
lines, in particular for the red [SII] lines. 

Using eq.[3]for replacing ij in eq.Q] we then have: 



< n, >= 



ki/k 2 - Gn 
G\iln\ -k 1 /(k 2 n 2 )' 



(4) 



This equation is next applied to the case of a stellar jet (see Fig. 
Q]), first considering spectroscopic observations with a long slit 
along the jet axis, where no spatial information across the jet is 
available. The observed line emissivity will be given by the in- 
tegration of the volume emissivity within the jet volume, where 
dV = 2nrdrAz, being Az the size of the beam in the z direction. 
If the electron density is represented by a Gaussian radial 



-r 2 /^ 2 



profile, then n e = n$e 
defined by eq.|4]becomes: 



and the spectroscopic mean density 



< n, >- 



(«1 - «2)(«0 ' 



J - n\L\ + n 2 L 2 



n\L\ - n 2 L 2 



(5) 



where L 12 = ln(l + «oMi,2)/(l + n ex /n 12 ), with «o the electron 
density on the jet axis and n ex the density evaluated at the jet ra- 
dius. This is the same expression derived by ABS89 for the case 
of an exponentially decreasing density in an isothermal stellar 
atmosphere. 

On the other hand, the mean electron density across the jet is 
given by: 



Jn e dV 



n - n e 



jdV log(n /« ex )' 



(6) 



where the right hand side term is obtained as before assuming a 
Gaussian density profile and dV = 2nrdrAz. 

The ratio between mean spectroscopic density and mean 
density (i. e. the ratio between the "observed" and the "real" den- 
sity) is shown in Fig. 2 as a function of the central electron den- 
sity for the [SII] ratio. Similar curves were calculated by ABS89 
for line ratios relevant to stellar atmospheres. 

As can be seen from Fig. [2] the ratio between the spec- 
troscopic and the mean density increases for steeper stratifi- 
cations of « e across the beam of the jet. At a fixed limb-to- 
centre electron density ratio n ex /no, the < n e > /« e ratio is 
slowly decreasing. The spectroscopic mean density is always 
larger than the mean density because denser regions contribute 
more to the emission. In the low density regime {no, n ex <K nj) 
in particular, it is straightforward to show that eq. [3] reduces to 
< « e > /n e « (2/3)(a 2 + a + 1) /(a 2 - 1) In a where a = no/n ex , 
and < n e > /« e ~ 2/3 In a when no » n eK as in the upper curves 
of Fig. m 

The results shown in Fig.|2]are particularly relevant for stel- 
lar jets because an important parameter that is determined obser- 
vationally is the mass flux (and the momentum flux, that has the 
same behaviour, as described below), defined by dM ~ pvrdr, 
where p is the mass density (oc « H = n e /xu) and v(r) is the jet ve- 
locity. If velocity and ionization fraction have a "top-hat" profile 
in the r direction, it may easily be shown that: 



M 



< M > <n e > 



(7) 



where M is the mass flux and the < M > is the one determined 
using the spectroscopic mean density. Therefore, values of mass 
and momentum-flux represent just an upper limit to the real val- 
ues (see e. g. Cabrit 2002]!. 
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Fig. 2. Ratio between the spectroscopic n sp and average density 
n av as a function of the density no on the jet axis, for different 
jet cross sections, characterized by their limb-to-centre electron 
density ratio n ex /no. 

2.2. Convolution effects for a resolved jet cross section 

We now consider the case where information on the emission 
profile across the jet is available. The volume and integrated 
emissivities i(r) and I(x) are related to each other by the Abel 
transform: 



I(x) = 2 



f 



i{r)rdr 
Vr 2 — x 2 



(8) 



where R is the jet radius and x is the projected separation be- 
tween a point in the observed image and the projection of the jet 
axis on the plane of the sky. 

In the case of an homogeneous medium i = j'o and eq.[8]leads 
to I(x) = 2io V/? 2 - x 2 . To show the effect of an inhomogeneous 
medium on plasma diagnostics, in the following we present re- 
sults obtained numerically integrating eq. [8] for different « e , T, 
xu profiles. 

Assuming a profile for n e , T, x H , for each radial position 
the synthetic emissivities in the [SII]/U6716, 6731, [NII]i6583, 
and [OI]/16300 lines are calculated, solving the statistical equi- 
librium equations using a 5-level atom (with atomic parameters 
from Mendoza 1983). These spatially dependent emission coef- 
ficients are then used to numerically integrate eq. [8] Finally, the 
physical parameters are determined from the integrated emission 
line intensities using the BE method (see AppendixlAl. 

The results obtained assuming a Gaussian profile for n e , and 
constant T and xu cross sections are shown in Fig. [3] The volume 
and integrated emissivities of [SII]/16716 and [SII]/i6731 versus 
r/R are shown in the top panel, both normalised to their values 
at r — 0. The [SII] ratio is also shown in the central panel of Fig. 
[3] while the volume and integrated electron densities are shown 
in the lower panel. 

The [SII] ratio determined using the integrated emissivities 
is higher than the original ratio (continuous and dotted line of 
Fig. 3, respectively). The difference between « e and A^ e increases 
from at the external radii to approximately 20% toward the jet 
axis. 

Fig. [4] shows the ratio no/No between the original and con- 
volved values of the electron density on the jet axis versus the 
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Fig. 3. Effect of the convolution on inhomogeneous jets. The 
dotted lines represent integrated quantities, while the continu- 
ous lines represent volume quantities. The electron density has 
a Gaussian profile ranging between 10 3 crrT 3 on the axis and 
100 crrT 3 at the jet radius. The temperature and the ionisation 
fraction are assumed to have constant values of 10 4 K and 0.1, 
respectively. Top: [SITJ/16716, /16731 integrated and volume in- 
tensities. Centre: [SII]/16716//16731 integrated and volume ratio. 
Bottom: volume and integrated electron density. 
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Fig. 4. Ratio between convolved and volumetric densities for dif- 
ferent values of «o/«ex- 



convolved electron density A^o, for different values of n e x/«o- 
As is evident from the figure, the convolution effect increases 
when the observed density approaches the critical density (of 



no and No goes up to ~ 60% for the most stratified case (bottom 
panel). 

Fig. [6] shows the results obtained for different initial strati- 
fications in n e , T and jch, corresponding to radial profiles with 
the form ~ I0~( r / R r where B — 2, 1,0.5 for the top, centre and 
bottom panels of Fig.|5]and|6]respectively. 

It is clear from Fig.|6]that effects similar to the one observed 
for « e in Fig. [3] are present also for T and xu- Additionally, Fig. 
[5] and [6] show that convolution effects increase with the plasma 
stratification (as expected). 

3. Data reconstruction 

3.1. Reconstruction of the jet structure using a Gaussian fit 
to the observed intensities 

The volume emission coefficient projected on the sky plane is 
also weighted on the beam size AxAz. Therefore, the observed 
emissivity Kx,) is given by: 



/(*,■) = 2Az 



r*i+\n r-R 
I dx I 



i(r)rdr 
Vr 2 - x 2 



(9) 



where i = -N, 0, ...,N, there are 2N + 1 observational values 
(with N = 0,1,2, .. .), and £±(#+1/2) = R- All the relevant quanti- 
ties are assumed to have negligible variations in the z direction. 

A simple scaling between i(f) and I(x) may be obtained in 
the case of an inhomogeneous medium, using the linearity of the 
Abel transform of a Gaussian distribution. In fact, if the volume 



emissivity has the form i(r) = i^e 



eq.|8]leads to: 



I(x) = yfno-i(x)erf -win — - 



-(!)' 



where 



erf(jc) 



2 



e ' dt 



(10) 



(11) 



is the error function, and z'o and / ex are the volume emissivity at 
the jet centre and limb, respectively. 

The property of the error function: erf(x) * 1 when x > 1 
implies that 



I(x) » -\[ncri(x) 

if x < R and ;'o » / ex ■ 

Using I(x) as given by eq.[T2] eq.|9]becomes: 



I(xi) = no- j'o (erf 



Xi+l/2 



Xj-l/2 ' 



erf- 



Az, 



(12) 



(13) 



If only one value of / is observed (jc !+ i/2 = — -Xj-i/2 = R)> the 
integrated emissivity will be simply given by: 

/ = Ko- 2 k) (l - e- R2/<rl ) Az ~ m^ioAz. (14) 

Given eq.[14l the relation between the integrated ratio G and the 
corresponding volume ratio g, in the case of Gaussian profile 
emissivities, is given by: 



10 4 cm -3 ) for the [SII] ratio. For high density stratifications, G - 18 where 



the no/No ratio converges to a maximum ratio, which depends on 
the value of A^o- 

Fig EJshows the results obtained for different profiles for n e 
and top hat cross sections for T and x^. The differences between 



<o-i 
T l = \ — 

\ cr 2 



(15) 



If o"i + cr 2 , the integrated and volume ratios will be different. 
Therefore, the physical parameters determined from the inte- 
grated ratio will be different with respect to the one obtained 
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Fig. 5. The dotted lines represent convolved quantities, and the 
continuous lines represent volume quantities. The curves corre- 
spond to the integrated and volume densities, obtained using the 
following initial profiles: n ~ 10 ^, T = 10000 K, x H = 0.1 
where ft — 2, 1, 0.5 in the upper, central and lower panel respec- 
tively. The volume electron density is 1000 cm -3 on the axis and 
100 cirT 3 at the jet radius. 
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Fig. 6. Effect of the convolution on inhomogeneous jets. The dot- 
ted lines represent integrated quantities, and the continuous lines 
represent volume quantities. Thin lines correspond to the ionisa- 
tion fraction, and thick lines to the temperature. The curves are 
obtained using the following profiles: T e , x, n e ~ 10~^ R ^, where 
jS = 2, 1, 0.5 in the upper, central and lower panel, respectively. 



using the volume ratio. Furthermore, the simple scaling given by 
eq. [15] may be used to reconstruct the volume intensities given 
the integrated intensities and the values of o~\ and <Ji, and from 
there to determine the radial dependence of the physical param- 
eters. 




Fig. 7. [SH]yl6716//l6731 ratio (points) and fits (lines) with eq. 
[TBI corresponding to different temperatures. The fitting param- 
eters are the following: n = 1841.89, 2457.4, 2883.9 cnr 3 , 
fli = 0.678, 0.687, 0.695, a 2 = 2.295, 2.333, 2.336 for T = 
5000, 10000, 15000 K, respectively. 

We first show how to apply eq. Q3] to the [SII] ratio to de- 
termine n e , and later the results will be generalised to more line 
ratios, to determine all of the physical parameters. 

A very good approximation to the relation between the [SII] 
ratio and n e may be obtained from eq.|4]and is given by: 

»e + no 

g= ; , (16) 

where the parameters no, a\ and a2 are in general temperature 
dependent. This formula is exact for a two-level atom, and in 
general is just a fit to the g = g{n e ) curve at constant T and xh. 
Fig. [7] shows the [SII] ratio and the fit calculated using eq. [T6l 
The fitting parameters are reported in the Figure caption as a 
function of temperature. 

In eq.[16] n e an d g represent the volume electron density and 
ratio respectively. A similar relation will also link the integrated 
quantities N e and G: 

G = jVe + "° . (17) 
N e a2 + noa\ 

Using eq.[T5][l6]and[T7]it is possible to obtain a relation between 
n e and A^ e : 



■ fli) + noa\{r\ - 1) 



«e = no— — -. (18) 

N e a2{i -rf} + n (a2 - V a V 

Fig. [8] shows the ratio between the volumetric and the observed 
densities n e /N e (obtained from eq.[T8l as function of A^ e , for dif- 
ferent values of r\. For electron densities A^ e S 5000 cm -3 and 
A^ e ^ 100 cm -3 , the n e /N e ratio becomes increasingly larger. 
These density ranges correspond to the high and low-density 
regimes (see Fig. |7). Also, for densities around 10 3 cirT 3 , the 
differences may be as high as a factor of 2 or 3, depending on 
the value of tj. 

This procedure may naturally be generalised to include the 
calculation of T and x^: 

- for each line with spatial information across the jet, the value 
of cr is obtained with a Gaussian fit to the observed emission 
profile. 
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- The volumetric line ratios are reconstructed using eq. [15] as 
gi,2 = h/h = o"2/o"iGi,2. 

- The reconstructed ratios are used to determine the physical 
parameters using any of the available methods (e. g. the BE 
method, see AppendixlAt. 

If the jet is observed through a narrow slit placed along the 
main jet axis, eq.Q~3]has to be used instead of eq.[14] Therefore, 
rj will be defined in this case as: 



/ cr^ 2 erf (Ax/cn) 
77 = erf (Ax/ crzY 

where Ax is the beam size along the x-direction. 



(19) 



3.2. Reconstruction of the jet structure using tomographic 
techniques 

The Abel transform (eq. [8) has been largely studied and applied 
in science, and complex methods were developed to solve it (e. g. 
Craig & Brown 1986 ). The analytical solution of the Abel trans- 
form is also well known: 



i(r) 



__i r R 

X Jr 



dl dx 



dx y/ x 2 _ r 2 



(20) 



The determination of z'(r) represents an "ill-posed" problem in 
the sense that small noise in I(x) produces large errors in the 
determination of z'(r), due to the numerical derivative present in 
the integral. Any numerical method solving eq. [20] using noisy 
observed data is therefore necessarily unstable. 

A method that is less affected by numerical instabilities is ob- 
tained by fitting the data by a smooth function and using eq.l20lto 
invert the fitted curve (e. g. Simonneau 1993). In this context, we 
use the multi-Gaussian method developed by Bendinelli (119911 1. 
which has been mainly used for deprojection of galaxy surface 
brightness distribution (e. g. Bendinelli 119911 Emsellem et al. 
[19951 Bendinelli & Parmeggiani [T9951 Cappellari 120021 . This 
method basically consists of producing a multi-Gaussian expan- 
sion fit of the data, and inverting the Gaussian series. Actually, 
this approach becomes very simple once the fit has been ob- 
tained, due to the already mentioned linearity of the Abel trans- 
form when applied to a Gaussian distribution. This implies that 
for a fit to a data series with a sum of Gaussians (centred on the 
same point): 



(21) 



;=1 



the Abel transform may be inverted analytically (eq.|20ll to give: 



1 ai 



-r 2 Jo? 



(22) 



where r is the distance from the jet axis. 



3.3. Application to HH30 

We now apply the techniques derived above to the HH30 jet. 
HH30 is an ideal candidate for this analysis. It moves nearly on 
the sky plane, has a clear side-to-side symmetry in the region 
close to the central star, and the cooling region is resolved spa- 
tially with Hubble Space Telescope (HST) observations (e. g. 
Burrows et al. IT9961 Ray et al. 119961 Bacciotti et al. 1999, 
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Fig. 8. Ratio of volume and integrated electron densities as a 
function of the integrated electron density for different values 

of 7]. 

Table 1. Data fit with Gaussian curves 



Gaussian Number 


x 


<x 


1 


20.97757 


153.4 


2 


20.98055 


30.62 


3 


20.98054 


16.42 


4 


20.98061 


15.71 


5 


20.98061 


15.69 



HM07). The data used for this analysis were taken by HST us- 
ing the slit-less spectroscopy technique, and presented by HM07. 
Using this powerful technique, the observations by HM07 re- 
solve the cross section of the HH30 jet with ~ 10-20 pixels. 

One position along the jet axis is considered here, at an an- 
gular separation of 3.9" from the HH30 source, corresponding 
to a region close to the central star-disk system. 

First, to calculate the centre xo of the Gaussian fit, i. e. the 
projection of the jet axis on the plane of the sky, all the observed 
profiles are added. The obtained emissivity is then fitted with 
an increasing number of Gaussians. The values of xq and the 
variance <x obtained using up to 5 Gaussians are shown in Tab. 
[TJ Fig- H] shows the data and a fit with 3 Gaussians (top), and 
the residual (bottom). From Tab. [T]and Fig. |9]it is evident that a 
very good approximation to the data is obtained already with 3 
Gaussians. 

The results obtained are shown in Fig.fTOland Fig.[TT] Each 
emissivity profile is fitted with a multi-Gaussian fit, and using 
eq. [22]the intensities as a function of the radial position, and the 
corresponding ratios are reconstructed. 

Fig.fTOlshows the observed data, the multi-Gaussian fit (top) 
and the reconstructed data (bottom). The ratios calculated using 
the reconstructed data differ from the original, observed ratios 
(see for example the [SII] /16716 curves with respect to the [Nil] 
in the upper and lower panels of Fig. ITOb . 

The physical quantities are shown in Fig. [TT] In the upper 
panel n e is shown. The points correspond to the observed values, 
while the continuous and dotted lines correspond to the multi- 
Gaussian fit and the reconstructed electron density values, re- 
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Fig. 9. Upper panel: Observed line emissivities (points) and 
multi-Gaussian fit as a function of the position across the jet us- 
ing 3 Gaussian curves. Lower panel: Residuals for a different 
number of multi-Gaussian fits. 



3,000- 
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Fig. 10. Upper panel: Observed line emissivities (points) against 
multi-Gaussian fit as a function of the position across the jet. 
Lower panel: Reconstructed line emissivity. 



spectively. The central and lower panels show the temperature 
and the ionisation fraction. 

Some interesting features may be noted from Fig. QT| First 
at all, it is evident that HH30 does not have a top hat cross sec- 
tion. The reconstructed n e , T and x H cross sections are much 
steeper than the measured cross sections. The reconstructed elec- 
tron density and ionisation fraction cross sections have strong, 
on-axis peaks, and the temperature has an on-axis valley. The 
total densities («h = n e /xu) are larger by a factor of two and 
four on the jet axis with respect to the jet radius, for the inferred 
from observation and reconstructed case respectively. The fact 
that we have lower temperatures in the denser regions is qualita- 
tively consistent with the stronger cooling that one would expect 
in higher density regions. The fact that the ionisation fraction is 
higher in the denser regions, however, is an effect that escapes 
simple qualitative arguments, and should be compared directly 
with predictions of jet formation and propagation. 

It is important to note that the original data were not decon- 
volved with the PSF (with a width of m 2.5 pixels), and we may 
expect a much stepper profiles using deconvolved data. Also, the 
jet presents some degree of asymmetry (see e. g. the wings of the 
observed temperature in the central panel of Fig-ITTTl. All of these 
effects, together with a complete study of the two-dimensional 
structure of the HH30 jet are left for a second paper. 

It should be noted that the results shown in Fig. 1 1 are af- 
fected by the fact that the [SII] ratio is close to the high den- 
sity regime (see Fig. 8). Other line ratios (see the discussion by 
HM07) should be used to deconvolve the jet beam closed to the 
source. 

4. Summary and discussion 

In this paper we have studied the influence of inhomogeneities 
on the electron density, temperature and ionisation fraction val- 
ues derived from observed line ratios. Additionally, we have pre- 
sented possible methods to analyse the data obtained from high 
angular resolution imaging of stellar jets. 

Simple scalings are possible when a minimum amount of in- 
formation is available. This minimum information is one value 
of emission line intensities for each position along the jet, to- 



gether with a determination of the jet width cr for all of the dif- 
ferent lines used for the plasma diagnostics. 

Bacciotti et al. (2000) observed DG Tau with 7 narrow slits 
parallel to the main jet axis, and at different positions across the 
jet, while Coffey et al. (|2007) used one slit placed perpendicular 
to the main jet axis. In both cases, they obtained spatial and kine- 
matic information across the jet axis. The approach presented in 
§3.1 can be easily applied to these data, calculating the average 
physical parameters and jet velocity as a function of position 
along the jet. The determination of the jet velocity has not been 
explicitly discussed in the present paper, but similar methods to 
the ones that we have discussed above can be used in this con- 
text. 

We have also shown that a more complex tomographic re- 
construction (see §3.2) can be used to reconstruct the cross sec- 
tion of the flow, obtaining a complete description of the three 
dimensional structure of the jet. 

The main limitations for the application of these techniques 
is that the jet axis has to lie close to the plane of the sky. 
However, it would be possible (but more complex) to reconstruct 
the 3D structure of the flow for arbitrary orientations of the jet 
axis. Additionally, the jet has to be nearly axisymmetric. Also, 
the applicability of the present deconvulation method is limited 
to the nearly axysymmetric region of the jet close to the source. 
Farther away from the source, in fact, the interaction of the jet 
with the ambient medium produces asymmetries in the intensity 
cross-sections. 

Finally, while this technique is currently applicable to a few 
observed objects, the use of future high resolution instruments of 
the new generation of telescopes will make it possible to apply 
this technique to a larger sample of objects. 

In this paper, preliminary application of our method to recon- 
structing the cross section of the HH 30 jet has been presented. 
In a second paper, we will describe a much more complete appli- 
cation of our method using observations of HH 30 and of other 
HHjets. 



8 



F. De Colle et al.: Diagnostics of inhomogeneous stellar jets 




20 25 
Radial position [pixels] 



Fig. 11. Values inferred from observations (points), fitted values 
(with a multi-Gaussian, continuous line) and reconstructed val- 
ues (dotted line) of electron density, temperature and ionisation 
fraction for the upper, central and lower panels respectively. 



Appendix A: BE method 

The emission coefficient (units: erg s _! cirT 3 sr _1 ) obtained by 
the transition from the i to j levels is given by: 



1 

An 



tij = —hvj iAi i 



rij n Ai n A 

ur-u "h, 

n Ai n A «h 



(A.l) 



where hvfj is the transition energy, Ay is the spontaneous 
Einstein coefficient, iij/n A ., n A ./n A and n A /n H are the excitation 
fraction, the ionisation fraction and the population fraction of the 
considered species. 

The BE method (from BE99) uses a series of line 
ratios to determine n e , T, x H : [SII]/16731/[SII]^6716, 
[OI]/16300/[SII]i6716 + 6731, [NII]^6583/[SII]i6716 + 6731. 
The [SII] ratio depends on n e and T, while the [OI]/[SII] and 
[NII]/[SII] ratios depend also on noi/nsn and nNii/nsn respec- 
tively. Additionally, nitrogen and oxygen ionisation fraction 
are determined assuming charge exchange equilibrium with 
hydrogen. Therefore noi/nsn and nNii/nsn become a function 
of the ionisation fraction and no/ns and n^/ns respectively. 
Finally, the sulphur is supposed to be all single ionised (n$n = 
ns) because the photoionisation rate due to diffuse UV radiation 
(e.g. Tielens|2005 pp. 267) is much larger than the SII radiative 
recombination rate. 

To calculate the physical quantities, we choose arbitrary val- 
ues of temperature and ionisation fraction, and invert the [SII] 
ratio deducing the electron density. Furthermore we use the N/S 
ratio and the electron density derived previously to determine 
the ionisation fraction, and finally we use the O/S ratio to find 
the electron temperature. Finally, we iterate until convergence 
using the new values of T and x^- 

The atomic parameters used to calculate collisional ion- 
isation, radiative + dielectronic recombination and charge 
exchange coefficients are from Cox ( 1 19701 ), Aldrovandi & 
Pequignot ( [19731 [19761 ) and Osterbrock ( fT989b . 
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